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ABSTRACT 

A linear stability analysis of ionized discs with a temperature gradient and an external 
axial magnetic field is presented. It is shown that both hydromagnctic and thermo- 
magnetic effects can lead to the amplification of waves and make discs unstable. The 
conditions under which the instabilities grow are found and the characteristic growth 
rate is calculated. The regimes at which both the thermomagnetic and magnetorota- 
tional instabilities can operate are discussed. 
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1 INTRODUCTION 



The origin of turbulence in astrophysical discs is often attributed to hydrodynamic and hydromagnetic inst abilities that can 
occ ur in differential l y rota ting stratified gas. The magn etorotational instability (MRI), first investigated by Velikhov ( 19591) 
and Chand rasckhar (jl960h and later fully recognized bv lBalbus fc Hawlevl (|l99ll) . is usually considered as one of the possible 



candidates to generate such turbulence and is thought to play an important role in the evolution and dynamics of astrophysical 
accretion discs. The growth of the instability in weakly ionized magnetized discs is of interest for models of star formation and 
the subsequent evolution of pro t ostellar discs. Numerical simulations of the MRI in accretion discs Jrlawlev. Gammie fc Balbusl 
19951 ; iBrandenburg et al.|[l995l ; iMatsumoto fc Taiimalll995l; iTorkelsson et al.lll99rj ; Urlt fc RiidigerlboOll ) show that turbu- 
lence generated can enhance essentially the an gular momen tu m transp o rt. T he MRI ha s been studied in detail f o r bot h 



stellar and accretion disc condit i ons (see, e.g., Fricke Jl969h : ISafronovl |l969h : lAchesonl (|l978h : iBalbus fc Hawlevl (jl99lh 



Kaisig. Taiima fc Lovelace! (|1992[ ); IZhang. Diamond fc Vishniad (|19940 ). In fact, the MRI can occur only in a relatively weak 
magnetic field, but a sufficiently strong field can suppress the instability completely (jBalbus fc Hawlevlll99ll ). This is related 
to the fact that the MRI is basically a long wavelength phenomena in the sense that it is stabilized for wavelengths shorter 
than A cr ~ 27tca/0 where ca is the Alfven velocity and Q is the angular velocity. Defining /3 = <? s /c?a where c B is the sound 
speed, we can estimate that the critical wavelength A cr is longer than the half-thickness of the disc, H ~ c 3 /0 for B > 1. 
More rigorous calculat ions for high and low values of the radial wave nu mber k r may be found i n ICoppi fc Kevesl <|2003h . and 
Liverts fc Mondl (|2009h . respectively. In particular, it has been shown bv lLiverts fc Mondl (2009) that the number of unstable 
MRI modes is decreasing with B. For example, it was found that there exist only three unstable modes in the disc for B ~ 10 3 . 
It is the goal of the present paper to study the effect of the temperature gradient on stability of magnetized discs in the case 
of high beta. It will be shown that the range of unstable wavelengths is significantly widened due to the vertical temperature 
gradients and can extend to values shorter than the disc thickness even in the case /3 > 1. 

In discs, the centrifugal force almost is balanced by the gravitational f orce, w h ile th e vertical structure appro x imately 
is det e rmined by hydrostatic equilibrium. The asymptotic analysis [see, e.g. lRegev (1983), iKluzniak fc Kital (|2000h . lOgiiviel 
( 19971 ). IShtemler, Mond, fc Riidigerl <|2009l) ] reveals that the vertical temperature gradient is comparable to that of the density 
in the case o f a sma ll aspect ratio. The des t abilizi ng effect of a tem perature gradient o n Alfve n waves has b een studies first 
by iGurevichl \ 19631 ): iGurevich fc Gel 'mond (|l967l ) and later on by IPoIginov fc Urpinl (|l979l ): lUrpinl (jl98ll ) for plasma in 
hydrostatic equilibrium. It has been shown that th e therm omagnetic instability (TMI) can arise that transforms a fraction of 
the thermal energy into magnetic one. Also, ICoppil (|2008T) has shown that the vertical temperature gradient in discs combined 
with the radial rotation shear give rise to vertically localized ballooning instabilities, the growth rate of which increases with 
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the temperature gradient. Thus, for a super adiabatic temperature gradient that can be obtained due to a strong heating 
source around the equatorial plane, the growth rate of the ensuing ballooning modes is found to be comparable to that of the 
"cylindrical" MRIs. 

In this paper, we consider the stability of hot magnetized discs taking into account the temperature gradient. We will 
take a look at the MRI and TMI in various ranges of the wavelength and determine the domains in a parameter space where 
each of those instabilities is dominant. 



2 BASIC EQUATIONS 



The dynamics of a magnetized disc is governed by the momentum, continuity, and induction equations: 
-Vp+ J x B + P G, 

i 

dp 



,dv _ ^> _ 
p(-+,.V« 



dt 

OB 
dt 



+ V ■ ( P v) = 0, 
•V x E, 



where G is gravity. 

The current density that is given by Ampere's law: 

V x B — p,oJ, 

where p,o is magnetic constant and B satisfies a divergence-free condition 
V-B = 0, 



(1) 
(2) 

(3) 

(4) 
(5) 



These set of equations should be complemented by the generalized Ohm's law that expresses the electric field E in terms 
of the electric current J and gradients of the thermodynamic quantities. It should be noted that due to quasi-neutrality of 
plas ma the ind e pende nt thermodynamic variables are only the temperature and pressure. Then, the Ohm's law reads (see, 
e-g. : 



Braginskiil < |l963h > 
1 



E = —v x B 



VP e + V ■ 3 + A • VT, 



(6) 



where the last two terms represent the galvanomagnetic and thermomagnetic effects, correspondingly. The tensors fj and A 
depend on the magnetic field B and, as a result, the transport properties of plasma are anisotropic with substantially different 
properties along and across the magnetic field if the field is sufficiently strong. The effect of the magnetic field on the transport 
properties is characterized by the magn etization param eter a c — ljbt where uib — qB/m c is the gyrofrequency of the electrons 
and r is their relaxation time (see, e.g.. lSpitzer ( 1978r)~). Ev en poorly conducting protostellar discs can be strongly magnetized 



if the electrons are the main charge carriers (jWardl 



1999h . 



The tensor terms in Eq.(6) can be written as follows 

0= r l J + r ll (Jx B) + r) 2 B(B- J), 
AVT = AVT + Ai(VT x B) + A 2 B(B ■ VT), 



(7) 

(8) 



where the second terms in both expressions represent the H all and the Nern s t effects, correspondingly. T he coefficients of 
these expressions can be obtained from the relations given bv lBraginskiil (|1963l ): iLifshitz fc Pitaevskiil (119811 ). It is convenient 
to represent the effects of the temperature gradients by introducing two quantities that have the dimension of velocity: 



Uit = AiVT 
u 2T = A 2 (B-VT) 



U 2 T 



(9) 
(10) 

ksT /m e \VT\ for moderate values of the magnetization parameter a e with ks 



both of which are of the order of uit 
being the Boltzmann constant. 

The basic state of the considered rotating discs is assumed to be characterized by an angular velocity that depends 
on the radial coordinate alone, O(r), and by hydrostatic equilibrium along the rotation axis z\ (r;ip;z) are the cylindrical 
coordinates. Thus, assuming polytropic equation of state, asymptotic analysis of the stead y-state properties of the discs 
demonstrates that the vertical tempera ture gradients are much larger than the radial ones (|Regev|[l98a : iKluzniak fc Kital 
200d : IShtemler. Mond. fc Riidigerl [20091) . Consequently, as a model problem, in the current calculations the steady-state 
temperature is assumed to vary only vertically. In addition, the background magnetic field is assumed to be a constant 
parallel to rotation axis z. 
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As a first step toward studying the TMI it is instructive to estimate the relative importance of the thermomagnetic 
and galvanomagnetic terms contribution to Ohm's law. For moderate values of the magnetization parameter and wavelength 
of perturbations measured by ca/O. (typical axial wave number for the MRI spectrum), assuming the disc thickness as a 
characteristic length scale for the temperature gradient one finds that the thermomagnetic terms in the Ohm's law and the 
induction equation are much bigger than dissipative ones arising due to resistivity terms if 

k B T 



2 2 „ 
UJ„T » 



(11) 



where lo p — yj n e q 2 /som e is the plasma frequency with n e being the number density of electrons and eo being the electric 
constant. This condition can be sat isfied for a very hot plasma with a rbitrary values of magnetic field and small or moderate 
number density. Thus for example, iBoettcher. Liang, fc Smith! (|1998f ) have investigated Galactic black hole candidates. Ap- 



plying their model to GX-339-4 the relevant data is that the disc temperature is of the order of 10 K with number density of 
the order of 10 24 m -3 , and a corona with temperature of the ord er of 10 8 K. Such range of parameters, with somewhat higher 
temperatures, has also been obtained by I Artemova et al.l (|2006r) . According to (JTTJ such conditions are indeed favorable for 
the onset of the TMI. The latter is consequently not suppressed by resistivity so that the galvanomagnetic terms can be 
dropped from Ohm's law. 

Even though condition (|11|) may be fulfilled in optically thick as well as in optically thin environments (as was exemplified 
in the previous paragraph) we restrict the current discussion to media such that radiative heat exchange is efficient in that 
medium, thus any thermal transport due to the electron conduction may be neglected and the energy equation may be 
dropped out. This indeed simplifies the subsequent calculations. To this end one may refer to optically thick medium where 
the radiative transport can be described as diffusive heat transport with a thermal diffusivity that depends on the temperature 
and therefore a temperature gradient can be well exist. It should be emphasized though, that the TMI may also be excited 
in optically thin discs, however in such environments, the energy equation should be taken into account. 

Consider a linear stability of the discs under axisymmetric short wavelength perturbations that propagate in ^-direction. 
The perturbed electric field E, magnetic field b, current density j and hydrodynamic velocity v vary in time and space 
according to exp(ifcz — iut) so local approximation is used. Retaining the thermomagnetic terms in Ohm's law ©, invoking 
incompressibility that results in zero pressure perturbations (the latter also means that the perturbations are restricted to 
Alfven modes) and assuming Keplerian angular velocity 0,(r) one can reduce linearized Eqs.(JT])-((6]) to the following system of 
equations: 



+ 3fi 2 - 



— U>a 



UJUJ2T ~ 



2iQ,U)2T —UJLU2T + 2iQuJiT 



-UJCJiT 



; iQ,U)2 



0. 



(12) 



Here the characteristic frequencies are the Alfven uia = kcA and the thermomagnetic oji^t = fcf i.2T frequencies. 

Before turning to the stability analysis it is instructive to notice that for long wave lengths perturbations all the matrix 
elements of second matrix in square brackets of the last system of equations are small and consequently the thermomagnetic 
effects can be neglected. In that case the MRI is the dominant mode of instability. We will elaborate on that point in the next 
section. The main interest of the current work however, is the case for which the matrix elements of second matrix in square 
brackets of the equations H12[) is comparable to the first one. This happens when: wi,2T ~ O, that indeed could occur in hot 
discs. 



3 DISPERSION EQUATION FOR THERMOMAGNETIC WAVES 

The dispersion relation that is obtained from Eqs.(|12p reads 



cj 4 + a-juj 3 + a20J 2 + a±uj + ao = 0, 



(13) 



where 



a3 = 

a2 = 
ai 

•i = 



2un 
~tl 2 



2w A + 



i — Q,U)2T 

2 



-2wiT(f2 +Ua) 



— 3£1 2 uja + t^A ~ 0, 2 (co 2 T + u)2t) + 2^ c * j 2t(3S1 2 — 5lj'a). 



(14) 
(15) 
(16) 
(17) 



If O = 0, Ea. (1131) reduces to the dispersion relation for Alfven waves mod ified by the temperature gradients (see Gurevicbl 

( 19631 ): iGurevich ^el'montl jl967h : iDolginov fc Urpinl (Il979h ; lUrpinl Jl98lh). If O ^ but VT = 0, Eq.^ reduces to the 



dispersion relation for the MRI obtained by 



Velikhov 



and 



Balbus fc Hawlevl (| 1 99 lh . which describes the combined 



effect of inertial and Alfven waves propagating in rotating discs. 
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Figure 1. The growth rate, obtained by numerical solutions of eg . (I13B for ut/ca = 0-2- 



It is instructive to examine a regime when the both MRI and TMI may operate. This implies the limit of thick disc 
(/3 S> 1) and uji^t 0- For typical wave number of the MRI spectrum (Q/ca) such regime reveals itself if £ = «2t/ca ^ 1- 
Furthermore one should note that in hot discs with the values of temperatures and number density quoted above the inequality 

can be hold for such regime and thus magnetic diffusivity may be neglected. A solution of the dispersion equation for 
£ = U2t/ca = 0.2 is depicted in Fig. 1. It is seen from this figure that there are two modes that become unstable. The 
upper curve corresponds to the MRI modified by the temperature gradient whereas the lower one describes a new instability 
that occurs due to the temperature gradient. It can be also obtained from the analytical consideration that properties of the 
instability are qualitatively different in the limiting cases of loa ^ fi (relatively long wavelengths) and to a f2 (relatively 
short wavelengths). In the long wavelength limit (oja ^ fi), the instability is predominantly the MRI with small corrections 
caused by the temperature gradient. In particular, the temperature gradient is responsible for a finite phase velocity which is 
vanishing in the isothermal case. In this limit, the eigenvalue is given by 

w = ±(5— +iV3au). (18) 
4 ca 

The corresponding eigenvectors are 

These eigenvectors represent the pair of linearly polarized waves that propagate in the opposite directions along the z-axis. 
The wave propagating in the positive direction (this direction has been assumed conventionally by setting a positive value for 
Q) is exponentially growing but the other one is decaying. 

The another branch of modes in the limit of relatively long wavelengths is characterized by oj ~ SI and small growth rate 
that is given by 

r = ImwM±|war^ (20) 
The eigenvalue and eigenvector for that branch are 

w = ±(n + il<* X ri£), b = {l, T 2i} T b, v = {l, l } T ^. (21) 
4 il A 2 J hop 



These eigenvectors describe the pair of left elliptically polarized waves that propagate in the opposite directions. In the absence 
of temperature gradients both waves are stable and represent inertial waves of the pure hydrodynamic nature. 

In the limit of relatively short wavelengths (at uj ~ uja fi) and for large temperature gradients such as ut/ca ^ 1 
the thermomagnetic effects play an important role in the development of instability. The Alfven waves can be significantly 
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modified by the temperature gradient and become unstable beyond the stability limit of the classical MRIs. The growth rate 
in this regime is given by 

r = c^^Z( 1± 2 ^ + " ) (22 ) 
2 Aloa 

Finally, in the short wa velength limit and under the c ondition ut /ca ^ 1 the rotation can be neglected and the problem 
reduces to that described bv lGurevich fc Gel'montl (| 19671 ). 



In order to gain deeper understanding into the nature of the TMI in the short wavelength limit, the dispersion equation 
is rewritten in terms of 6 = b x ± ib y in the limit — > 0. Then, the dispersion equation reads 

[lj 2 + u){u)\t =F iw2r) — Ua\ b ± = (23) 

This equation describes four circularly polarized waves. One pair of waves propagates along the magnetic field, and its 
frequency is 

u = u A 2~ ~2~' ^ ' 

The corresponding eigenvectors are 

6 ± = 6, & T = 



Wit i . W2T 

— 1 — ± i- 

2oja "2oja 



j t = 0. (25) 

i>0 



CA 

The unstable wave (upper sign) is right circularly polarized while the damped wave (lower sign) is left circularly polarized. 

The wave for which the phase shift between perturbations of the magnetic field and velocity is smaller than tv (it is the 
appropriate value for pure Alfven waves) grows exponentially. Another wave for which the phase shift between the magnetic 
and velocity perturbations is larger than n is damped. The other pair also comprises of two right and left polarized waves 
that propagate in the opposite direction to the magnetic field. One of these waves grows exponentially whereas the other is 
damped in the same manner as in the first pair. Summarizing, right circularly polarized waves propagating in the opposite 
directions along the magnetic field grow exponentially with the growth rate 

r = (26) 

The picture of circularly polarized unstable waves allows an additional interpretation of the instability. To that end it 
should be noted that j r± = ^kb* /(j,q where j* is the complex conjugate amplitude of the current density of the circularly 
polarized waves. The thermomagnetic contribution to the perturbed electric field is anti Hermitian as it can be seen from 
Eq. (|25[) . Then, it is easy to show that the average work over the wave period is given by 

{E.j)=-^\b\ 2 . (27) 
The time derivative of the average energy (W) for the wave propagating along the external magnetic field is given by 

This equation together with Eq. <|27p clearly shows that the wave energy grows exponentially with the rate 2I\ 



4 SUMMARY 



We have considered the thermomagnetic instability in ionized non-isothermal Keplerian discs that are threaded by an ex- 
ternal magnetic field. It has been shown that, in the short wavelength limit, the temperature gradient allows to develop the 
thermomagnetorotational instability for wavelengths satisfying the condition cja ^> for which the MRI does not occurs. 
The growth rate of this instability can be comparable to that of the MRI. Physically, the thermomagnetorotational instability 
transforms a fraction of the th ermal energy car ri ed out by the heat flux into the magnetic energy of Alfven waves (modi fied by 
the temperature gradient, see iGurevichl (|1963h : iGurevich fc Gel'montl (|1967h : iDolginov fc Urpinl (|l979h : IUrpinl jl98lh ). The 
presence of shear, however, changes substantially the growth rate of instability. The considered combined magnetorotational 
and thermomagnetic instability can play an important role in dynamics of astrophysical discs. 
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